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Measurements of the energy spectrum and of the vortex-density fluctuation spectrum in superfluid 
turbulence seem to contradict each other. Using a numerical model, we show that at each instance 
of time the total vortex line density can be decomposed into two parts: one formed by metastable 
bundles of coherent vortices, and one in which the vortices are randomly oriented. We show that the 
former is responsible for the observed Kolmogorov energy spectrum, and the latter for the spectrum 
of the vortex line density fluctuations. 
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Below a critical temperature, liquid helium becomes a 
two-fluid system in which an inviscid superfluid compo- 
nent coexists with a viscous normal fluid component. The 
flow of the superfluid is irrotational: superfluid vorticity 
is confined to vortex lines of atomic thickness around 
which the circulation takes a fixed value n (the quantum 
of circulation) . Superfiuid turbulence [U |2] is easily cre- 
ated by stirring either helium isotope (''He or ■^He-B), 
and consists of a tangle of reconnecting vortex filaments 
which interact with each other and with the viscous nor- 
mal fluid (which may be laminar or turbulent). The most 
important observable quantity is the vortex line density 
L (vortex length per unit volume), from which one infers 
the average distance between vortex lines, i ~ L~^l^ . 
Our interest is in the properties of superfluid turbulence 
and their similarities with ordinary turbulence. 

Experiments [3l|4] have revealed that, if the superfluid 
turbulence is driven by grids or propellers, the distribu- 
tion of the turbulent kinetic energy over length scales 
larger than I obeys the celebrated k~^/^ Kolmogorov 
scaling which is observed in ordinary (classical) turbu- 
lence. Here k is the magnitude of the three-dimensional 
wavenumber (wavenumber and frequency are related by 
k — f /v, where v is the mean flow). Numerical calcu- 
lations performed using either the vortex filament model 
[3 [3] or the Gross-Pitaevskii equation [71 [S] confirm the 
Kolmogorov scaling. Another important experimental 
observation is that in both ''He 9J and '^He-B |10], the 
frequency spectrum of the fiuctuations of L has a de- 
creasing /~^/^ scaling typical of passive objects [HI [TT] 
advected by a turbulent flow. This latter result seems 
to contradict the interpretation of L as a measure of su- 
perfluid vorticity, u) = kL which is usually made in the 
literature [2 [5]. In fact, from dimensional analysis, the 
vorticity spectrum corresponding to the Kolmogorov law 
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should increase with / (as /^Z"^), not decrease. Since the 
vortex line density is a positive quantity, a better anal- 
ogy is to the enstrophy spectrum: however in classical 
turbulence this spectrum is essentially flat [TH [13], in 
disagreement with the helium experiments [HI 110] . 

The aim of this letter is to reconcile these two sets of 
observations (each separately backed by numerical simu- 
lations). We shall show that, at any instant, the vortex 
tangle can be decomposed into two parts: vortex lines 
which are locally polarised in the same direction, form- 
ing metastable coherent bundles, and vortex lines which 
are randomly oriented in space. The former is responsi- 
ble for the Kolmogorov energy spectrum, and the latter 
for frequency spectrum of the vortex line density. 

Following Schwarz [TJ], we model vortex filaments as 
space curves s(^, i) which move according to 

ds 

— = + as' X (v„ - Vs) - a's' x (s' x (v„ - v,)) , (1) 
at 

where t is time, a and a' are known temperature depen- 
dent friction coefficients [T^, s' = ds/d^ is the unit tan- 
gent vector at the point s, ^ is arc length, and v„ is the 
normal fluid velocity at the point s. We set the temper- 
ature to T = 1.9 K, typical of many finite temperature 
studies (corresponding to a = 0.206 and a' = 0.0083). 
The self-induced velocity of the vortex line at the point 
s is given by the Biot-Savart law [IF 

where k = 9.97 x 10^'' cm^/s (in "^He) and the line in- 
tegral extends over the entire vortex configuration C. 
The calculation is performed in a periodic cube of size 
-D = 0.1 cm. The numerical techniques to discretize the 
vortex lines into a number of points (j ~ 1,---N) 
held at minimum separation compute the time 

evolution, de-singularize the Biot-Savart integrals, eval- 
uate Vs using a tree-method (with critical opening angle 
0.4), and algorithmically perform vortex reconnections 
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FIG. 1. (Color online) The time evolution of the total vortex 
line density (black dashed line) L (cm~^) and the polarised 
component (cm~^) (solid red line). 

when vortex lines come sufficiently close to each other, 
are all described in our previous papers [HI [T7] . Here we 
take = 2 x 10"'^ cm and a timestep of 5 x 10~^ s. 
The number of discretization points varies with time; 
in the statistically steady state which we shall discuss, 
N Ki 4000. 

The turbulent normal fluid is modelled by the following 
synthetic turbulent flow |18) : 

M 

v„(s, i) = ^(A„xkmCOS(/)m + BmXk„sin0m), (3) 

m— 1 

where 0,„ = k,„ • s + /,„i, k„ and /,„ = y/k^E{km) 
are wave vectors and angular frequencies. This flow is 
solenoidal and time-dependent; with a suitable random 
choice of A^, and k^ (adapted to the periodic box). 
The normal fluid's energy spectrum has Kolmogorov 
form E(k„i) ^ km in the range from ki (corresponding 
to the integral scale) to fcjv/ (corresponding to the dissipa- 
tion scale). This synthetic model of turbulence compares 
very well with Lagrangian statistics obtained in direct 
numerical simulations of the Navier-Stokes equations and 
experiments. Here we take M = 188, ki = 0.1 cm and 
kM = 1.8 X 10~^ cm, which corresponds to the Reynolds 
number Re„ = {ki/kM)'^^^ ~ 200. For computational 
simplicity our model ignores any back-reaction of the vor- 
tex lines onto the normal fluid. 

We integrate the vortex lines in time according to 
Eq. ([T]), for a period of 10 s (approximately 25 large eddy 
turnover times of the normal fluid). We find that, after 
an initial transient, the vortex line density saturates to 
a statistically steady state, as shown in Fig. [l] indepen- 
dently of the initial condition (various vortex loop con- 
figurations were tried). 

We analyse the properties of the superfluid turbulence 
in the statistically steady state {t > 1.2s). Firstly, we 
compute the frequency spectrum of the fluctuations of 
the vortex line density about its average value < L >= 
1.15 X 10** cm~^. Fig. [2] shows that the spectrum scales 
as f~^^^ for large /, as observed in experiments [HI HO] 
and numerical simulations [B] . Secondly, we compute the 
superfluid energy spectrum E(k), defined by 





FIG. 2. (Color online). Power spectral density, PSD, (ar- 
bitrary units) of fluctuations of the total vortex line density 
L (black upper solid line) and of the polarised vortex line 
density (red lower solid line) vs frequency / (s~^). The 
dashed (magenta) line shows the f~^^^ scaling. 

where V = is volume. The top curve (labelled a) 
of Fig. |4] shows that the energy spectrum (computed on 
a 256'^ Cartesian mesh) is consistent with the cascad- 
ing classical Kolmogorov scaling E{k) ~ k~^^^ in the 
range kn = "in/ D < k < kg = 2Tr/i, in agreement with 
experiments [31 [3] and numerical simulations [5H5]; at 
larger wavenumbers k > kf the spectrum exhibits the 
non-cascading scaling of individual vortex lines. We con- 
clude that the numerical model thus reproduces the two 
major spectral features (the energy spectrum and the 
vortex density fluctuations spectrum) observed in super- 
fluid turbulence. 

To solve the puzzle described in the introduction, we 
examine the homogeneity of the turbulence. A simple 
three-dimensional plot of the vortex lines may give the 
wrong impression that the vortex tangle is spatially uni- 
form: the orientation of the lines is not apparent and 
some lines partially hide each other. A more careful anal- 
ysis is required. We deflne a smoothed vorticity field ijJ 
at the discretization points using a kernel with finite 
support, the M4 kernel [TS] (effectively a cubic spline): 

N 

u;(s,) = «^s;.iy(r,;„/i)AO, (5) 

where = |s, - Sj|, A^^ = |sj+i - Sj|, W{r,h) = 
g{r / h) / {-Kh?) , h is a, characteristic length scale, and 

{l-lq^ + lq\ 0<q<l; 
g{q)=l\{2~q)\ l<q<2; (6) 

[0, q>2. 

This approach is commonly employed in the smoothed 
particle hydrodynamics (SPH) literature [19j . and means 
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FIG. 3. (Color online) Left: Snapshot of the vortex tangle £ 
local magnitude of the smoothed vorticity field uj. Right: ' 
u){si) > 1.4a;rms) and random components (bottom, u){si) < 
maximum value cjmax = 38.4 s~^; cjrms ~ 15.7 s~^ 

that we only need to take the contribution of discretiza- 
tion points within a radius of 2h from each point s^. We 
tested this smoothing algorithm against a Gaussian ker- 
nel g{q) = e~'^ [3D] and found comparable results for 
single snapshots of the vortex tangle. The advantage of 
the new algorithm is that u) is computed only at the dis- 
cretization points, so it can be evaluated during the time 
evolution. 

Setting h — £, the smoothed vorticity field lj al- 
lows us to identify the presence of any coherent vor- 
tex structures. From the vorticity values u}j at each 
discretization point s^- we compute the rms vorticity 

Wrms = \/^^^il7l^^^^Ksj)P- and, following Ref. [11], we 
decompose the vortex lines into locally 'polarised' and 
'unpolarised' fields, i|| and Lx respectively. The po- 
larised vortex line density consists of the discretization 
points Sj associated to intense metastable vortical regions 
where the magnitude of the smoothed vorticity exceeds 
the rms vorticity by a threshold value, uJ^Si) > lAuj^^s 
(we shall discuss this threshold shortly). The remaining 
discretization points form the 'unpolarised' field Lx- As 
with the total vortex line density, we find that Ly rapidly 
saturates to a statistically steady state, see Fig. [ij 

A snapshot of the vortex tangle in which the vortex 
lines are coloured according to the local value of the 
smoothed vorticity is shown in Fig. |3] (left); vortex bun- 
dles are clearly visible. The right-hand side of Fig. [3] 
shows respectively the polarised (top, uj{sj) > lAuj^i^s) 
and unpolarised vortex lines (bottom, uj{sj) < 1.4a;i-ms): 
the former consists of bundles and contributes to Ly, the 
latter is spatially random and contributes to Lx- 

We find that the polarised and unpolarised vortex lines 
have different properties. The frequency spectrum of the 
fluctuations of Lx scales as f~^^^, as for the total line 
density L (top black curve of Fig. [2]), whereas the spec- 
trum of the fluctuations of L|| slowly increases at low fre- 



: t = 7s. Vortex lines are locally coloured according to the 
'he same snapshot but split into the locally polarized (top, 
1.4a;rms) respectively. The color scale is normalised by the 

quency and drops rapidly at high frequency (bottom red 
curve of Fig. [2]). Changing the threshold value within 
the range 1.4 to 2.2uj^jns does not change this distinc- 
tion. The flattening and the slow increase with / of the 
spectrum of the fluctuations of the polarised vortex lines 
are consistent with what is observed in ordinary turbu- 
lence [12, 13 . 

We now turn our attention to the role of the polarised 
lines on the energy spectrum. We have seen that the top 
curve (labelled a) of the top panel of Fig. |4] displays, in 
the range kjj < k < ki, the Kolmogorov energy spectrum 
corresponding to all vortex lines. The curves labelled b, 
c, d, and e show energy spectra arising only from dis- 
cretization points Sj such that w is below a given (de- 
creasing) rms threshold. It is clear that, as we remove 
the contribution of the high-intensity bundles, the energy 
spectrum E{k) in the range kjj < k < kg becomes shal- 
lower, and changes from a k~^^^ scaling to a fc~^ scaling. 
The bottom panel of Fig. [4] compares the energy spec- 
trum arising from the isotropic vorticity Lx (top curve, 
labelled a) and that arising from the polarised vorticity 
Lj| (bottom curve, labelled b). It is apparent that the vor- 
tex bundles correspond to the classical k~^/^ Kolmogorov 
scaling, and that the random vorticity corresponds to the 
k~^ spectrum. 

The decomposition of L into L|| and Lx is robust and 
holds during the time evolution, thus vindicating Roche 
and Barenghi who suggested it when discussing an 
experiment [S]. Unlike that experiment, in our calcula- 
tion the random vortex lines contain most of the energy. 
A possible explanation of this difference is that our model 
of turbulent normal fluid does not contain strong vorti- 
cal structures, unlike real turbulence or DNS, thus our 
polarised bundles are an underestimate of reality (the 
presence of vortical structures in the normal fluid would 
certainly induce stronger vortex bundles in the super- 
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FIG. 4. Energy spectra E{k) (arbitrary units) vs wavenum- 
ber k (cm~^) corresponding to Fig. [3] Top: The upper solid 
line (a) is the energy spectrum corresponds to the flow induced 
by all vortex lines. Lower curves correspond to the flow in- 
duced by vortices with smoothed vorticity below the following 
thresholds: (b) uj < l.Tojrms, (c) lAuJrms, (d) 1.2LJrm.s, and 
(e) LUrms- Bottom: Energy spectra corresponding to vortex 
lines with smoothed vorticity respectively below (a) and above 
(b) the threshold lAuJrms- The dashed lines display the 
and the k~^^^ scalings. 



fluid, as shown in numerical simulations [5T], hence in- 
crease the energy contained in Ly). It is also interesting 
to remark that the the smoothed vorticity u) can be be 
interpreted as the coarse-grained vorticity of the HVBK 
equations j22|. and that the decomposition of L into 
and Lx, has an analogy to that of Farge et al. [23] for 
classical turbulence. 



In summary, we have shown compelling evidence that, 
at any instant, the turbulent vortex tangle can be de- 
composed into a polarised and a random component. 
The polarised component, associated to iy, consists of 
metastable coherent vortex bundles, shown in the top 
right panel of Fig. [3] and is responsible for the observed 
Kolmogorov k~^/^ energy spectrum. The random com- 
ponent, shown in the bottom right panel of Fig. |3] is 
responsible for the observed /^^^^ frequency spectrum 
of the fluctuations of the vortex length. The result not 
only solves an apparent puzzle between experiments, but 
justifies theoretical attempts such as Lipniacki's [24] to 
develop a two-scale approach to superfluid hydrodynam- 
ics. 
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